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In this work, we consider the Casimir effect due to massless fermionic fields in the presence of 
long cylinders. More precisely, we consider the interaction between a cylinder parallel to a plate, 
between two parallel cylinders outside each other, and between a cylinder lying parallelly inside 
another cylinder. We derive the explicit formulas for the Casimir interaction energies and compute 
the leading and the next-to-leading order terms of the small separation asymptotic expansions. As 
expected, the leading order terms coincide with the proximity force approximations. We compare 
the results of the next-to-leading order terms of different quantum fields, and show that our results 
support the ansatz of derivative expansions. 
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I. INTRODUCTION 

Casimir effect is one of the most interesting predictions of quantum field theory that has been verified experimentally. 
In [l|, Casimir proposed that there exists a force between two parallel perfectly conducting plates due to the vacuum 
fluctuations of the electromagnetic field confined between the plates. The idea of Casimir energy is natural. In 
quantum theory, the ground state energy of a simple harmonic oscillator is not zero, but is equal to &u/2, where w is 
the angular frequency of the associated oscillator. A quantum field can be considered as the superposition of infinitely 
many simple harmonic oscillators with different frequencies. This led Casimir to define the Casimir energy to be 

the sum of the ground state energies of the quantum field. This is an infinite sum that needs regularization. However, 
without the presence of boundaries or external conditions, the Casimir effect would not be manifested. The Casimir 
effect is most interesting in the presence of two objects, such as parallel plates. In principle, after subtracting the self 
energies, there should be some finite amount of energy left, which would create interaction between the two objects. 

The first experiment that successfully verified the presence of Casimir effect appeared near the end of the 20th 
century [^. This has stimulated another surge in the research activities in Casimir effect, especially in conjunction 
with the development of nanotechnology. Under the same reasoning, the definition of Casimir energy (HD works not 
only for electromagnetic fields as originally proposed by Casimir, but for any quantum fields. However, for fermionic 
fields, one should add a minus sign in front of the formula due to the different spin and statistical behavior of fermionic 
fields. In fact, such definitions of Casimir energies have been used to compute the Casimir effect of two parallel plates, 
and the Casimir self energies of spheres, cylinders, etc, in the second half of the 20th century. From the point of view 
of statistical physics, such a definition is natural since the Casimir energy so defined appears as the zero temperature 
part of the free energy. 

Even though there is intensive research in Casimir effect since 1980s, for a long time, it is not clear how to 
compute the Casimir interaction between two objects, except by using approximations. Around 2006, various groups 
of researchers simultaneously tackled this problem for some particular geometries using quantum field theory methods 
such as Green’s functions, path integrals, wave expansions, etc, that can be more or less categorized as multiple 
scattering approach UHli or mode summation approach [TthI^ . The general method for arbitrary objects has been 
synthesized in [2^ for scalar fields and i n pH for electromagnetic fields. Both of these papers approach the problem 
using multiple scattering formalism. In [22l| , we used mode summation approach to interpret the formulas derived in 
^ 1^ . An advantage of our formalism is that it is not restricted to scalar fields or electromagnetic fields, and it is 
also not restricted to (3-|-l)-dimensional Minkowski spacetime. For example, we have used the formalism in to 
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compute the sphere-plate and sphere-sphere Casimir interaction in {D -|- l)-dimensional Minkowski spacetime in 
and [^ . 

The studies of Casimir effect of fermionic fields can be dated back to 1980s. Before the end of the 20th century, 
there are a few works that considered Casimir effect of fermionic fields [2^ - 1^ . but this is a relatively small number 
compared to the research works in Casimir effect of scalar fields and electromagnetic fields. One of the possible reasons 
is that ferminoic fields are relatively harder to deal with. Nevetheless, after entering the 21st century, Casimir effect 
of fermionic fields has started to attract attention and there are a number of works in this area [29l - l4q . with some 
applications to carbon nanotubes. 

In ( 4 ^ . we considered the fermionic Casimir interaction between two spheres and derived the small separation 
asymptotic behaviors. For application to carbon nanotubes, it is natural to consider the cylindrical geometries. In 
this work, we consider the fermionic Casimir effect between a cylinder and a plate, and between two parallel cylinders. 
For parallel cylinders, we consider both possible cases: one is where the two cylinders are outside each other, and 
one is where one cylinder is inside the other. We derive the explicit formulas for the Casimir interaction energies and 
compute the small separation asymptotic behaviors. The results are compared to the results of other quantum fields. 
We also use our results to stipulate the ansatz of derivative expansions proposed in [s^. 


II. THE CASIMIR INTERACTION ENERGY 


A. Plane waves and cylindrical waves 


In this work, we consider the Casimir interaction between a cylinder and a plate, and between two cylinders due to 
the vacuum fluctuations of a massless Dirac field ip which satisfies the equation 


= 0 . 


( 2 ) 


Here = 9/^ + F^, and F^ is the spin connection. 

On the boundaries of the cylinders or plate, we impose the MIT bag boundary conditions: 


(1 -h i^^n^)ip 


boundary 


= 0. 


(3) 


To derive the Casimir interaction energy, we use the formalism we developed in [22j| . 

First we need to solve the equation of motion ([3]) in rectangular and cylindrical coordinates. 

For cylinders, we align them so that their axes of symmetry are parallel to the z-direction. Then using the cylindrical 
coordinates 


X = pcosip, y = pshup, z = z, 


(4) 


a cylinder of radius R can be described as p = i? if its axis of symmetry is the z-axis. 

When considering the cylinder-plate interaction, we will assume that the plate is given by a: = L, where L > R, so 
that the cylinder is parallel to the plate. The plane waves are then parametrized by (ky, kz). The fermionic waves 
can be divided into positive energy modes and negative energy modes, as well as regular waves and outgoing waves, 
each has two families: They can be written as 


—2sgn* kxX-\-ikyy-\-ikzZ^iijjt 


where 




4 (±).* 

ky kz ,2 


v± 

/ 


T 


I 

0 

-sgn^kx +iky 

k 

0 
1 

sgn^kj; +iky 


i 
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(5) 


(6) 
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Here k = —, * = reg or out, sgn,,gg = 1, sgu^^^ = — 1 and 


U — _ 1.2 _ U2 

— V ^ ' 

13 5 

In cylindrical coordinates, the fermionic waves are parametrized by m and kz, where m = ± —, ± —, ± —,..., 


„/,(±)>* _ r* r(=*=)"* JkzZ^ibjt 




n(±).* _ k V 


r(±)i* _ k\ ^ Y l^ 

^r„fc..2- T^r_,{k^p)e<^-^> 




A:_l = x/kP' — kl, 


r:^{z) = uz), = 


^reg _ ^out _ 

^ 5 '-'m 2 ^ 

are normalization constants introduced to facilitate the change to imaginary frequencies. 


B. The Casimir interaction energy between a cylinder and a plate 


Assume that the cylinder has radius i?, length H and its axis of symmetry is the z-axis. Let the plate be described 
by a; = L, with center at O' = (L, 0,0) and dimensions H x H. Let x' = x — L, L = Lex- 

In the region between the cylinder and the plate, the Dirac fields can be represented in two ways: one is in terms 
of the cylindrical coordinate system centered at O: 


m— — oo \ j = l,2 


/OO ’—' / _ _ 

-00 ™ \ — 1 O • — T o 


and one is in terms of the rectangular coordinate system centered at O': 


=H> r^ cf>-‘*'“43 ;"(x',„)]. (n) 

J —00 j —00 j —00 \ X_1 O X_1 o / 


These two representations are related by translation matrices V and W: 


d5;7(Y.u.)\ ^ /yyy \ ^sy(x.„) 


kyk ^,2 V 
j (±),out 




Ci:2‘(x,.a)J 7-00 27r ’ 
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so that 


/{±),mk^ 

^^(ih) ,fcy kz '' 

jl±),kykz 


=H 


°° tik 

dKy / *rrL,k. ^ \ / Ci 


= E 


/_^ 27r y , 

\ m^ky m,ky 

^ ky ,m ky.m / 


(±),21 -,r (±),22 I I (±),kykz I > 


'^{±),mkz 

,{±),mkz 


m— — oo 


The boundary conditions on the cylinder give a relation of the form 

/i{±),mkz\ ( {^),mkz 

I E -t(±) (“i 

1 ^(±),m/c2 I mkz \^{^),mkz 
In imaginary frequency, uj = k = in, ^ = \/) 

TLt=±* 





where 


si = kR {K^_^{ciR)I^_^{nR) - if^+i (7-R)Wi(7-R)) ■ 


The boundary conditions on the plate gives a relation of the form 


^ {:t),kykz\ / j(=l=)?^yfez 

l±),kykz I kykz I ,l±),kykz 


where 


fw -^1 

^ h h -T 

ri-yKz ^ 


\jl'^ + kl + ky ik^ 

-ikz Jl’^ + k^-ky 


For the translation matrices V and W defined by m, using techniques introduced in we find that 

/I 


' 1 ky 


0 




As discussed in 


where 


TT 

n 


n _ o_ 

\/T+E '+k 
1“ - 





l^ + kl 


the Casimir interaction energy is then given by 


ii(±) ^j{±) jj 

m.m' mkz 


2yj. m,ky kykz ky^rr 


, I - . m-\-m' 

+ kl + ky\^ g-2^7=+fe?i 


K \-ikz -7 


7 —jfcj 


dky 

27r 


The integral can be computed explicitly to give 

/ = =F -Tmfc 

m,m' ~ _ ^rriKz 

r\j 


7 —ikz 

-ikz -7 


R^m+m' ( 27 A) ■ 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

( 20 ) 

( 21 ) 

( 22 ) 

(23) 
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C. The Casimir interaction energy between two cylinders 

Consider two cylinders of length H and radii Ra and Rb respectively. The axis of symmetry of the cylinders are 
given respectively by a; = y = 0 and x = L,y = 0, both parallel to the z axis. 

We consider two scenarios: 

I. The two cylinders are outside each other. In this case, L > Ra + Rb and d = L — Ra — i?_B is the distance 
between the two cylinders. 

II. The cylinder of radius Ra is inside the cylinder of radius Rb- In this case, L < Rb — Ra and d = Rb — Ra — L 
is the distance between the two cylinders. 

In the region between the two cylinders, the Dirac fields can be represented in two ways: one is in terms of the 
cylindrical coordinate system centered at O: 


=iJ / dw ^ ^ 

m— — oc \_7 — 1,2 


{±),mk^j^(±),Teg 


At7(='.“)+ E if 

1 = 1.2 


(24) 


and one is in terms of the cylindrical coordinate system centered at O': 


/ oo OQ / _ 

E E 

m' — — oo yj — 1,2 


/ (±),reg/ / 




(25) 


i=i,2 


The two representations are related by translation matrices. In case that the two cylinders are outside each other. 


vsird-'.")), f- 


fez ,1 

),OUl 

fez,2 

In case one cylinder is inside the other, 

/,(±),i-eg/ / 

/(±),OUt/ N 

(x,a;) 


rr(=l=)i2lN 
^m,m' ^ m,m' 

m— — oo m,m' ^m,m' / 


^fj{±)M ^C±),21^ 
^ m'm ^ m'm 


= E 




lie I I L I I L .lit 

fr(±)A2 jj{±),22 

i ^ m ^ m m ' m J 


m— — i 
oo 

= E 


''I'mkz ,1 

(±),out/ N 
mkz ,2 (X,W) 




^y(±).ll 

m,m m,m \ 

y{±)A2 y{±),22 I 

m— — (X) ^m,m' / 

m' ,m m' ,m 

m' = —oo \ m',m m',m / 




vV’^^2S®(x,w)y ’ 

U^fe°2‘(x',w)y 


(26) 


(27) 


Using the method in [^, one finds that 


UL-I,, =(-i)^ 


CL 

uS™ =(-!)' 


/(±) 


=(-l)^ 


^ ra^ —m{f 
^ mi —m{f 
-^Im-m'{lL)l, 


(28) 


As in the case of electromagnetic fields, we find that the translation matrices are all equal to a scalar times the identity 
matrix. 

When the two cylinders are outside each other, solving the boundary conditions on the two cylinders give the 
matrix as in m- The Casimir interaction energy is given by 


= Ell>"(•-““’«>) 


(29) 
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with 


=TLt (Ra) 

m" 


( 30 ) 


When the cylinder with radius Ra is inside the cylinder with radius Rb, the boundary condition on the cylinder 
with radius Ra still give the same as given by dig. However, the boundary conditions on the cylinder with 
radius Rb gives 


{^),m'kz 


= -T^ , 

m' 


(±) ( 4 




kz ( I ’ 


(31) 


where 


Tpmfe. =±i. 


-iT' 

ikz 
7 


ikz 
7 
1 =F 


>^Rb [iI_i{iRb) + II^^{iRb)) y 

with given by dig. The Casimir interaction energy is then given by the same expression ( 1291 ) but with 

{Ra) Y. {RB)^^5'',m' 

m" 

=TL1 {Ra) E {RB)Im"-m'{lL). 


(32) 


(33) 


III. SMALL SEPARATION ASYMPTOTIC BEHAVIOR 

Casimir effect will be most significant when the separation between the objects is small. However, the computation 
of the small separation asymptotics of the Casimir interaction energy is often a tedious problem. Nonetheless, a 
systematic method has been developed in a series of papers Eliil, [ 23 , [ 23 , IHlI - lb^ . Making the substitution 

kz = u cos a, K = It sin a 


dsg and dig, we find that 


RCas — 


hcH 

47r^ 


udu f daTr Eln (I-MW) . 

^ 0 _j_ _ 


Expanding the logarithm and trace, we have 
hcH 1 


47r^ ^ ^ s + 1 

+ .-s=0 


udu / da E] • • ■ E/ 


momi mim2 rrismo 


rriQ mi iris 


Here the trace tr is the trace over 2x2 matrices. 

In the following, we are going to discuss the different cases separately. 


(34) 


(35) 


A. The cylinder-plate case 


uj = Ru, 




In the cylinder-plate case, let 
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Then 


[(i) 

rrii 


1 

ujsiT? a(Kl^_,{u:) + 


-1 =F iS§i sin a i cos a 
i cos a 1 =F sin a 


1 —icosa 

— l cos Of —1 


K, 


rrii+mi+i 


(2w(l + £)) 


— sin^ a =F i^i sin a sin a cos a 

sin a cos a — sin^ a ± sin a 


Kmi+mi+i (2w(l + e)). 


where 


Let 




(36) 


uiQ = m, rrii = m rii, 


m 


UJ = 


VT^' 


Then the main contribution to the Casimir interaction energy comes from m~£ ^, rii ^ e ). 
Using the Debye asymptotic expansions of modified Bessel functions: 


I^{vz) 


K^{vz) 


1 


i^ri{z 




\/2t:v (1 + z'^)i 


i 1 + 


TT e 


-vri{z) 


2j^(1 + z2)4 


1 - 


V 

uiit{z)) 


where 


(37) 


ri{z) =\/l + z'^ + log 
t{z) = 


1 + VTT^' 


vTT^’ 


”■<*>=8^ M' 


(38) 


one can check that the term has order e, and hence would not contribute to the small separation asymptotic 
expansion up to the next-to-leading order term. Hence, 


[(i) 


1 




K, 




[2oj{l + £))) I. 


(39) 


This is independent of a and -I-, —. Moreover, ignoring terms with order higher than e, 
Hence, after replacing summation by corresponding integrations. 


— 

rrii+i — 


^ hcH ^ 1 dr 


dmm^ f dm... f dn,tr J 


(40) 


Now 


Krrii+mi+i (2u;(l -|- £)) 




W 1 -f 


mi-\-4 


(w)' 


1 (w) 

m.A — 4- k / 


I. 


Using the Debye asymptotic expansions, one can find the asymptotic expansions for 

Kmi+rUi+i (2a;(l -|- £)) 

m.z — i ^ 2 


(41) 



and 


1 + 


mi-\- 


l{uj) 


1 (w) 

m.A — 4- ^ '' 


separately, up to terms of order e. These give an expansion of the form 





T{ni 


4m 


(1 + -^1,1 + ^ 4 , 2 ) I, 


(42) 


where Ai_i and Ai ^2 are respectively terms of order y/e and e. Ai^i is an odd function in and and hence 
integrating it over an even function of rii gives 0. As a result, the next-to-leading order term in the small separation 
asymptotic expansion of the Casimir interaction energy is of order e smaller than the leading order term. We have 


2hcH ^ (-1)"+^ 1 

7ri?2 S -I- 1 

8=0 


3 — 5 

dTT 2 


dmm 2 


2'*+r7r~ Jo 


dni... dus 


X exp — - 


2e(s-|-l)m T(ni - 


(43) 


-E 


4 = 0 


4m 


2^0 j—2+1 


2 = 0 


The integration over rii are standard Gaussian integrations and the integration over m and r can also be performed 
explicitly. We find that 


Eca 


ihcH 


16V2^£5i?2 (s + 1) 


E 


(- 1 ) 


s+1 




36 


R 


I'K^hcH'/R 


3840V2d 

One can easily check that the leading term 


T_ 

'm 


20 


2l7r2 


3840V2dt 


(44) 


(45) 


coincides with the proximity force approximation. 

The correspon ding results for Dirichlet (D), Neumann (N) and perfectly conducting (C) boundary conditions have 
been obtained in [l0l |: 


E 


D 

Cas 


E 


N 

Cas 


E 


c 

Cas 


TT^hcEly/R f 

1920y2di V 

TT^hcE[\/R / 

1920y2di V 

TT^hcEly/R f 

960V2di V 


7 d\ 


36 


■ 7 

40 


37r^ 

■ 7 

20 





(46) 


B. The case of two cylinders outside each other 


In the case that two parallel cylinders of radii Ra and Rb are outside each other, let 


uj = {Ra + Rb)u, 
Ra 

a =- 

Ra + Rb 


e = 


Ra + Rb 
Rb 


h = 


Ra + Rb 
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Then 


where 


Let 


1 ^m'^+rrii ((1 “t” ((1 T 

wsina(iL^__i(aa;)+iL^^^i(aw)) 6wsina M + M) 

— 1 =F sin a i cos a \ f— 1 ± sin a i cos a 


i cos a 1 =F i£§i sin aj \ i cos a 1 ± i'^i sin a 


=auj (^K^_i{auj)I^_i{auj) - K^.^i{auj)I^.+ i{au})j , 
% =buj (^K^,_i{buj)I^,_i{buj) - i (&a;)/„,+ i (few)) . 

mo = TO, mi = m + rii, 

,b bb 

mi = ■^{rrii + mi+i) + qi = + rii+i) + Qi, 


UJ = 


ly/l — t'^ 


(47) 


(48) 


As before, to has order e rii and qi has order e =. 

Now, 

— l=F*^iSina icosa \ /—lii'^jsina icosa '\ _ /sin^ a ± i — %) sin a ± (.^i — “^i) sinacosa 

icosa l=Fi^jSinay icosa lii'^jsinay ± (^i — "^i) sinacosa sin^ a =F * (-^i — “^0 sin a 

(49)' 

In the same way, we find that the ^i and ‘^i terms would not contribute to the leading and next-to-leading order 
terms of the Casimir interaction energy. 

Hence, 


fc -| p\ 7 p OO n OO p OO 

^Cas TZPY 73 J dni... J dn.ti (m^ ^ J , 


ttR^ —- 

^ s—0 


where 


^rrii ,mi+i 


((1 “1“ ((1 H“ „ 

(iQi - v li. 




(A'-b-U"*') + AJ„+. (m^)) 

Using the Debye asymptotic expansions (|37p again, we can find the small e expansions for 

((1 -|- e)uj)Kmi^ 1 +mJ ((i + g)^) 

1 (aw)Ar^ , 1 (buj) ’ 

m,- —-4 ^ ■' m — 4 '' '' 


(50) 


(51) 




1-f 


K^, i{buj) 

m.' — i A ' 




1-f 


1 (aw) 

m..- — i A ' 


and 
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up to terms of order e. These give 


■^mi ,mi+i 


27rm J_ 


f 2em brim - a^r 2'\ mt 

dqi exp ( ——--—g, ) (1 + Bi^i + Bi^2) li¬ 


ar 


4m 


After the Gaussian integration over qi, we have 


^rrii ,7712 4-1 


'Jin 

2-jTTm 


exp - 


2£m bririi — rii+i)" 


ar 


4m 


(1 + Cip + Ci^2) I- 


Hence, 


Eca 


2hcH g (-1)"+^ b^ 


s-3 


n<X) 


■kR^a 


s=0 


s + 1 2®+i7r 2 Jo 


drr 2 / dmm 2 


dni... 


dn„ 


X exp — - 


2e(s+l)m bT{ni — nj+iY 


J —00 

S —1 S 


ar 


2=0 


4m 


1+E E 


2=0 


As before, integrations over rii, m and r give 

^ ihcHai'Jb ^ (-1)"+^ (, 1 

^ W'j2^eiR\ (^TTF V ■ 


s=0 

Ttt^ HcH\JRaRb 


1 - 


12 Ra + -Rs 
7 d 


'7 (s + 1)^'' 
36 9 


■ + 


7 20 


The leading term 


3840V-Ra + Rb V 12 Ra + Rb 

hcH\JR aRb 


36 2l7r2 


d d 
Ra Rb 
d d 
Ra Rb 


Eca 


3840V2diVR-4 +Rb 


(52) 


(53) 


(54) 


(55) 


(56) 


also coincides with the proximity force approximation. 

The corresponding results for Dirichlet, Neumann and perfectly conducting boundary conditions were obtained in 
[ 5 ^. They are given by 


EEa. 


rpN 

^Cas 


Ecas 


1 - 


n^hcH-^RARB 
1920V2dtVR-4 + RB 
TT^hcH^/RARB 
1920V2dVR a + Rb V" 12 Ra + Rb 

d 


12 Ra + Rb 36 \Ra Rb 




1 - 


7 


■k^HcEI'JRaRb 


1 - 


7 


■ + 


■ 7 40 

7 20 


d d 
Ra Rb 
d d 
Ra Rb 


(57) 


960^2^5 VRa + Rb V 12 Ra + Rb 

Notice that for all different boundary conditions, the ratio of the next-to-leading order term to the leading order term 
contains the universal terms 

7 d 
12 Ra -I- Rb 


and 


7_ _d_\ 

36 \Ra ^ Rb) ■ 
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C. The case of one cylinder inside another 

In the case that the cylinder of radius Ra is inside the cylinder of radius i?_B, let 

w = {Rb - Ra)u = 


TO\/l — 


£ = 


QT 

d _ ^ L 
Rb ~ Ra R’ 


Ra 


Rb — Ra ’ 


b = 


Rb 


Rb — Ra ’ 


, b , , b & - 

=m + ni, = —(n^ + Ui+i) + qi = -m+ —(rii + ni+i) + Qi- 
la a la 


Then as in the case of two cylinders outside each other, we find that 

1 


rrii ,mi+i 


, ^m'—rrii ((1 ^m'-—rrii^i ((1 £')^) 

dqi —^^^-1. 


J-co bu (^Il,_,{bu;) + 


In the same way, we have 

^ mcHaiVb ^ (- 1)^+1 f, 

'^“^16V27r£ii?2 (s + l)4 Iv 


+ 


7 


'-A s=0 

Ttt^ hcH \/RaRb 
3840V2di^RB - Ra 


1 + 


12 Rb — Ra 

7 d 
12 Rb — Ra 


7 (s + l)2i 
36 9 


+ 


7 


20 


36 2l7r2 


_d _^ 

Ra Rb 

jL a. 

Ra Rb 


(58) 


(59) 


As computed in the corresponding results for Dirichlet, Neumann and perfectly conducting boundary conditions 
are given by 



TT'^hcH^J RaRb 


Cas 

imoA^diy/RB- Ra 


TT’^hcH^/ RaRb 

A 

Cas 

l920V2di ^/Rb - Ra 

'C 

'n'^hcHy/ RaRb / 

'i+- 

Cas 

960^2^1 V-Rb - Ra ' 


7 


7 


12 Rb — Ra 36 \Ra Rb 


7 


12 Rb — Ra 
7 d 


12 Rb — Ra 


+ 


7 40 

L 


d d 
Ra Rb 
d d 
Ra Rb 


(60) 


Again, we notice that for all different boundary conditions, the ratio of the next-to-leading order term to the leading 
order term contains the universal terms 


and 


12 Rb — Ra 


7 

36 \Ra Rb 


IV. DISCUSSIONS 

In the previous section, we have seen that the leading order terms of the Casimir interaction energies are indeed 
equal to those derived using proximity force approximation. The asymptotic expansion of the Casimir interaction 
energy up to the next-to-leading order term has been a subject of much interest. In [59l - [^ . Tosco et al performed 
derivative expansion to the path integral representation of the Casimir energy and obtained an integral expression for 
the expansion of the Casimir energy up to the next-to-leading order which is not completely convergent. However, 
it is good enough for computing the next-to-leading order term. Their method works successfully for scalar field but 
only partially for electromagnetic field. No results have been derived so far for fermionic fields. 
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Inspired by the work 
expansion of the form 


Bimonte et al proposed in [50l| that the Casimir interaction energy has a derivative 




• Vi/i + 

J 2 

+ /3x (-ff)Vili • Vila + P- {H)z ■ (Viii X Vila) + ..•), 


(61) 


where is the Casimir energy density between two parallel plates, E can be taken to be the z = 0 plane parametrized 
by X = {x, y), z = i?i(x) and z = i?a(x) are the height profiles of the two objects with respect to E, and H = Hi — H 2 
is the height difference. The leading term 




is precisely the proximity force approximation. Using the invariance of the Casimir interaction energy with respect 
to tilting the reference plane E, it was found that 

P-{H) =0, 

/3x(-ff) =2-/3i(iJ)-/3a(i?). 

For Dirichlet, Neumann and perfectly conducting boundary conditions, /3 = /3i = /3a is found to be a pure number 
that only depends on the boundary conditions, which is given by 


3^=-n- 


Since 


= 3 ''- 

2 
■3 


£LiH) = 


30 

TT^ 

15 

-TT^ 


a 


iJ3’ 


where a depends only on boundary conditions, if one let 


(62) 


(63) 


Co = / 


cii = / d^x- 


VHi ■ VHi 




f .2 VFla • VFla 
f .2 veil • VFla 


then the ansatz m says that 

^Cas ^ Q^|co + P{cil + C 22 — 2ci 2) + 2ci2|. 
For cylinder-plate configuration, one finds that 


Co 


_3ttVRH / d 

~ 4:V2di V 4i? ■ 


?)TT'/RH 2 d 

~ 4V2d^ ^ 3)R 
C22 =Cl2 = 0. 


(64) 


(65) 


( 66 ) 
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Therefore, the ansatz dm) says that 


pDE 

-^Cas 


37r/Ri/ r, /2^ 1\ d 


( 67 ) 


As observed in the results from exact computation (H51) do indeed satisfy (1571) with the various /3 given in (1571) . 
Our new result (l44l) satisfies (I67|) if the value of /3 for fermionic fields with MIT bag boundary conditions is given by 


For two cylinders outside each other, 

SttH \/RaRb 


Co = 


Cll = 


C22 = 


Cl2 = 


4:\/2d^^^ VRa + Rb 
3ttH \J Ra Rb 


/3 = o- 


1 + 


10 


4:\/2d^/‘^ \/Ra + Rb 

3t:H \J RaRb 

4\/2d5/2 7^^ + Rb 

3ttH \J RaRb 


4 Ra + Rb 

2 d 

3 Ra T Rb 
2 d 


4\/2d®/^ V Ra Rb \ 3 Ra + Rb 


3 Ra T Rb 
2 d 


1 
4 

2 d 

3 Ra 

2 d 

3 Rb 


jL 

Ra 


Hence, (I55|) gives 


E. 


DE 

Cas 


3'!t^/RaI^H 

AV2d^VRA+RB' 


a < 1 — 


7 


12 Ra + Rb 




d 

Rb 


_d_ 

Ra 


d 

Rb 


( 68 ) 


(69) 


(70) 


One can check immediately that our new result ((55|) satisfies (1701) with fip given by 

A similar computation for two spheres also shows that our results in [i^ satisfies (105|) with Pf given by (I68p . 
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